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Strain qualitatively changes the low-energy band structure of bilayer graphene, leading to the 
appearance of a pair of low-energy Dirac cones near each corner of the Brillouin zone, and a 
Lifshitz transition (a saddle point in the dispersion relation) at an energy proportional to the 
strain [M. Mucha-Kruczynski, I.L. Aleiner, and V.I. Fal'ko, Phys. Rev. B 84, 041404 (2011)]. 
Here, we show that in the vicinity of the Lifshitz transition the conductance of a ballistic n-p and 
n-p-n junction exhibits an anomaly: a non-monotonic temperature and chemical potential depen- 
dence, with the size depending on the crystallographic orientation of the principal axis of the strain 
tensor. This effect is characteristic for junctions between regions of different polarity (n-p and n-p-n 
junctions), while there is no anomaly in junctions between regions of the same polarity (n-n' and 
n-n'-n junctions). 

PACS numbers: 73.22.Pr, 62.20.-x, 71.70.Di 



I. INTRODUCTION 

Bilayer graphene (BLG)j^2l a crystal consisting of 
two graphene monolayers arranged according to Bernal 
stacking^ is a material with versatile properties. In con- 
trast to monolayer graphene, where the linear dispersion 
(Dirac cones) near each corner of the Brillouin zone (K 
and K' points) is very difficult to alter, the low-energy 
band structure of BLG can be qualitatively modified by 
relatively weak external perturbations. For example, a 
transverse electric field opens a mini-gap in the BLG 
spectrumpJ^Sl Also, it has been shown that a relatively 
small uniaxial strain (of only a few percent) leads to 
a change in the topology of the low-energy dispersion, 
which then exhibits two Dirac mini-cones near each cor- 
ner of the Brillouin zone (see Fig. llJPfflD Both in the 
conduction and valence bands, these cones are connected 
by a saddle point at which the Fermi lines reconnect, 
a configuration which is known as a Lifshitz transition 
(LiTr)pJII2] (l n contrast, in a monolayer homogeneous 
strain only results in a small shift of the Dirac cones 
away from the corners of the Brillouin zone, without any 
qualitative change of the linear dispersion or the chiral 
properties of the electrons.') 

In this paper, we study transport characteristics of an 
ideally clean homogeneously strained BLG crystal, aim- 
ing to find features in the temperature and chemical po- 
tential dependence of its two-terminal conductance that 
would reflect the presence of the saddle point in the dis- 
persion relation. We consider a short and infinitely wide 
strained BLG strip (the only geometry where strain in a 
two-terminal device would be homogeneous^), adjacent 
to BLG regions suspended over metallic contacts. The 
information encoded in the two-terminal conductance of 
such a device of finite length is complementary to what is 
manifested by the sheet conductivity of an infinite flake 
discussed in Ref. [TS]. Since contacts with metals heav- 
ily dope graphene, we model the BLG terminals with a 
high (e.g., n-type) density of carriers, whereas the strip 
in the middle is considered to be at a low density of car- 




FIG. 1. (Color online) Left: Top view of an unperturbed 
(top panel) and a strained (bottom panel) bilayer graphene 
(BLG) crystal. The top and bottom layers are shown in yellow 
and red, respectively. Strain modifies the intra-layer nearest 
neighbors coupling 70, as well as the inter-layer coupling 73 
between atoms at the center of the other layer's hexagons. 
Right: Electronic band structure in the vicinity of the Bril- 
louin zone corners K and K', with focus on the low-energy 
dispersion near the K point for unperturbed and strained 
BLG. 

riers (either of n- or p-type, which can be controlled by 
an external gate). We choose the amount of strain in 
the structure such that it induces a LiTr at the energy 
of about ±5meV, measured from the charge neutrality 
point. According to Ref. [5], such an effect on the bands 
can be generated by about ~ 1% of uniaxial strain. Note 
that in suspended graphene structures" strain of such 
size may be inflicted involuntarily, either by processing 
and annealing of the flake, or by displacements of con- 
tacts due to the different contractions upon cooling of 
the substrate and the supporting metallic electrodes. 

Our findings show that the dependence of the conduc- 
tance G(fi, T) on the chemical potential and tempera- 
ture does indeed reflect the spectral reconstruction by 
strain, in the form of a conductance anomaly: a non- 
monotonic dependence of G(/i, T) on both parameters, 
fi and T . This behavior is characteristic for the regime 



where the chemical potential is close to the saddle point 
in the electron/hole spectrum on one of the sides of an 
n-p junction, or in the middle of an n-p-n device. The 
conductance anomaly is sensitive to the crystallographic 
orientation of the sample geometry, which determines the 
principal axis of the strain tensor. These results are de- 
scribed in detail in Sections |III|a nd |lV[ with their graphic 
representation shown in Figs. [3] and [4] The calculations 
are based on the Landauer-Biittiker approach,^ with 
transmission probabilities obtained in the transfer ma- 
trix methodpSKH] Section [n] introduces the model for a 
strained bilayer graphene device and identifies the prop- 
agating and evanescent modes required for these calcula- 
tions. 



II. ELECTRON DISPERSION AND 
PROPAGATING MODES IN STRAINED BLG 

In this section we identify the energy dispersion and 
transport modes in homogenously strained BLG regions. 
These results are used in the subsequent sections to study 
the transport in devices made out of several such regions 
(n-p and n-n' junctions in Sec. Ill n-p-n and n-n'-n junc- 
tions in Sec. 



IV I 



The lattice structure and parametrization of the mini- 
mal rele vant tight-binding model for electrons in strained 
BLCP^ are illustrated in the left panel of Fig. [l] The 
stacked layers have every A site within each layer sur- 
rounded by three B sites and vice versa, with intralayer 
coupling 70 ~ 3eV; A 2 sites are on top of Bi sites, 
with interlayer coupling 71 ~ 0.4eV, while A1/B2 sites 
sit over/under the hexagons in the other layer and are 
coupled by 'skew' hopping energy 73 ~ 0.3eV. The low- 
energy electronic states reside on the sites Ai and B 2 , 
while the sites A 2 and B\ support states in split bands 
which do not contribute to low-energy transport. For 
unstrained BLG the low-energy states near each corner 
of the Brillouin zone form two approximately parabolic 
bands, a valence band and a conductance band, which 
touch each other at the K or K' point, as shown in the 
top right panel of Fig. [T] 

Uniaxial strain changes the intralayer and interlayer 
hopping integrals 70 and 73 by making them direction 
dependent, as shown in the bottom left panel of Fig. [T] 
Neglecting trigonal warping for large enough strain, the 
corresponding low-energy dispersion near a given cor- 
ner of the Brillouin zone is described by the effective 
HamiltoniarPiSl 



V(x) 



4K> +' 

V{x) 



-2iij> 



(1) 



Here m ~ 0.035m e is the effective mass, tt = p x + ip y 
parameterizes the in-plane momentum relative to the 
K or K' point, and we" 2 "* accounts for the change 
of the couplings due to the strain, where (f> is the an- 
gle between the principal axis of the strain tensor and 




FIG. 2. (Color online) Schematic representation of a sus- 
pended BLG device with strain axis oriented along a;-direction 
(as defined in Fig. [TJ. The sketch illustrates the example 
of an n-p-n configuration of such a device (fi < 0). In the 
highly doped contact regions the Fermi level (dotted line) lies 
high up in the conduction band (yellow), where the disper- 
sion is parabolic. In the central region the Fermi level lies in 
the valence band (red), and is close to the charge neutrality 
point, where the dispersion is modified due to the two Dirac 
mini-cones and the saddle point associated with the Lifshitz 
transition. Shading indicates occupied states. 



the crystallographic direction of the crystal. Using the 
tight binding model for BLG, one finds®^ that w = 
(3/4) (773 - ?7o)73(<5 - witn r /o,3 = din 70,3/d In r^s 
where tab is the distance between carbon sites, while 5 
and 5' are the two principal values of the strain tensor. 

Near each corner of the Brillouin zone, the low-energy 
dispersion relation obtained from Eq. (JlJ exhibits two 
Dirac mini-cones, which are separated from the parabolic 
spectrum at high energies (w <C |e| < 7i/2) by a saddle 
point at e = ±w (see bottom right panel of Fig. [I]). For 
energies |e| < w between the saddle points each mini-cone 
results in a disconnected, approximately circular Fermi 
line. At the saddle point the lines connect pair wise in a 
LiTr, and beyond the LiTr there is only a single Fermi 
line encircling the K or K' point. Relative to these corner 
points, the strain-induced Dirac points are positioned in 
the momentum plane at 



Po = Po ( cos 4>, sin < 



Po 



= ±V2r 



(2) 



After expanding H in Eq. ([TJ in momentum p— po around 
these Dirac points (and keeping only linear terms), we 
find that each is characterized by a Dirac velocity v* = 
Po/ '(2m). 

Below, we study how this strain-induced change in the 
topology of the electronic bands affects the transport 
properties of a device made of a strained flake of BLG, 
where a narrow and wide strip of width L y 3> L x sits be- 
tween two highly doped BLG regions suspended over two 
metallic contacts. The band alignment in such a device is 
sketched in Fig. [2j Metallic electrodes heavily dope BLG 
in the vicinity of the contacts, thus determining two leads 
(regions I and III) with a high carrier density. An exter- 
nal electrostatic gate controls the doping, and, thus, the 
chemical potential fj, of the electrons in the middle part of 
the flake (region II) , which we consider to be close to the 
neutrality point. In our model this doping profile is taken 
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into account by potential steps at the sample edges, 
if 0<i<L x , 71 



V(x) 



-Vo 



otherwise, 



> V > w, (i. 



In the remainder of this section, we identify the transport 
modes in the various regions of the system. 

The stated conditions make the energy dispersion in 
the leads approximately parabolic, e ps p 2 /(2m) — Vo, 
and the plane wave states the same as chiral states in 
unstrained BLG,^ with a very little effect of the strain. 
For given incidence angle 9 of an incoming electron in 
the contact, we parameterize their transverse momentum 
along the step as p y = y/2m\Vo + ij sin(0), and use the 
longitudinal component 



y l £ = I ^n2m\V a 



l = ±, 



(3) 



to characterize propagating modes (n = +, real momen- 
tum) and evanescent modes (n = — , complex momen- 
tum). For propagating modes in the conduction band of 
the leads, as considered here, the group velocity is di- 
rected parallel to the momentum, and thus the index I 
coincides with the propagation direction along the x axis, 
i.e., I = + denotes a state propagating to the right. 

We now turn to the modes in the weakly doped re- 
gion II. The left panels in Fig. [3] show isoenergetic lines 
for electrons in the valence band at low energies for un- 
perturbed (Fig. [3^ a)) as well as strained BLG (Fig.^b, 
c, d)), with w — 5meV and for several orientations of 
the strain principal axis. States corresponding to plane 
waves moving to the right are indicated by red and green, 
and to the left by purple and blue. (Note that for some 
of these modes the group velocity is directed opposite to 
their momentum.) These isoenergetic lines reflect that 
the low-energy dispersion relation in strained graphene 
is determined by the modified condition 



1 



(Pi + Py 



m 



4m 2 W« ' ry, m W [tf x ~ P 2 y ) COS(20) 

— 2wp x p y sin(20) + w 2 



(4) 



(here and in the following " ~ " denotes quantities specific 
for region II; energy e and transverse momentum p y are 
conserved for elastic scattering at a straight interface). 
For given values of energy and transverse momentum, 
this equation may have four, two, or no real solutions 
p x , where the latter situation arises at any given fixed 
energy beyond a critical value \p y \ — p y , c , which depends 
on the orientation of the applied strain. We denote the 
corresponding propagation direction in region I by 9 C , 
p y c — y^mjVo + ej sin(0 c ), which signifies the critical 
angle beyond which electrons from the lead only couple 
into evanescent modes. Below the critical value and for 
large values of |e|, there are always two real and two 
complex solutions, while for small |e| there are two or 
four real solutions, which depends on the propagation 
direction and on the orientation of the applied strain, as 
we now discuss in detail. 



For the unstrained case [neglecting w in Eq. @], the 
parameters in region II are given [in analogy to Eq. (|3|] 
by 



fx = I sign(e) 



^Jn2m\e\ — p 2 



(5) 



Here, p x + (p x + ) i s real, arL d corresponds to right- 
moving (left-moving) plane waves, while p~~ (p x ~~) cor ~ 
responds to evanescent waves decaying to the right (left). 
(The factor sign(e) accounts for the fact that in the va- 
lence band the group velocity is directed opposite to the 
momentum.) 

For strained BLG with strain orientation = 0, we 
find from Eq. Q that 



- I sign(e) 



P7 



Am 2 t 



8mwp 2 



pi + 2mw, (6) 



where n,l — ±. The left panel of Fig. (3^b) shows ex- 
amples of several isoenergetic lines, with strain-induced 
Dirac points on the axis p x in the momentum space. By 
inspecting Eq. ^ one notices that for |e| < w and p y < 
y/me 2 /(2w) [angles where sin(0) < y / e 2 /(4w|V r + if)] all 
four momenta are real [p x + (red),p~~ (green), p+~ (pur- 
ple), p~ + (blue)] and the Fermi line is split into two 
pockets. When e is slightly below (above) the LiTr in 
the valence (conduction) band, |e| > w, the Fermi line 
is continuous but deformed. For small values of \p y \, 
Eq. Q then gives two real solutions (pj + , P x + ) and two 
imaginary solutions {p x ~ , Px), while for larger values 
of \p y \ (just below the critical value p y , c ) there are four 
real solutions. 

Fig. [3jc) illustrates the propagating modes for strain 
with orientation <fi — 7r/4, where the momenta were found 
numerically from Eq. Q. The four colors correspond to 
the right-moving plane waves (p x + red, p~~ green) and 
left-moving plane waves (p x ~ purple, p~ + blue). We 
now find at most two real solutions for fixed energy and 
transverse momentum. Above the LiTr in the valence 
band, there is a range of transverse momenta around p y = 
(normal incidence from the leads) in which there are 
no propagating modes in region II. 

For the strain axis oriented at (f> = 7r/2 [Fig. J3jd)] , the 
four solutions of Eq. Q are 



p x = I sign(e) y n y 4m 2 e 2 + 8mwp 2 — p 2 — 2mw, (7) 

where, as before, n = ± and I = ±. By inspecting Eq. ^ 
we find that for all energies and angles below 9 C only the 
momenta p^ x + are real. In Fig. |3jd), the corresponding 
propagating waves are marked red (pj + , right-moving) 
and blue (p x + , left-moving). Above the LiTr in the va- 
lence band, there is again a range of transverse momenta 
around p y — in which no propagating modes exist in 
region II. 
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III. TRANSPORT ACROSS N-P AND N-N' 
JUNCTIONS 

In this section we study how the presence of strain 
affects the electron transmission across a single potential 
step, from a heavily doped region I to the low-density 
region II. Depending on the sign of the doping, this can 
be an n-p (/i < 0) or n-n' (a > 0) junction. In the 
Landauer-Buttiker approach paHU ^he conductance G of 
such a junction is determined by the energy and angular 
dependence of the transmission probability T(e, 0) of an 
electron in the conduction band incident from the left to 
emerge in the valence (or conduction) band at the right 
of the interface. 



A. Transmission probability 

To calculate the transmission probability T(e, 9) we 
employ the transfer matrix method.^! Using separation 
of variables (allowed for a straight interface) the spinor 
eigenstates of Eq. |l]) can be written as $i,n(x, y) = 
^1,11 (x) e lPyV , where 




and 

^ = v»Ve ' 

Pin = • 

e 

Here, as before, indices I and II label regions to the 
left and right from the potential step, I, n = ± dis- 
criminate the branches of longitudinal momentum, a; n , 
bi n are the wave amplitudes, and vi n = l/h\de/dp l ™\, 
Vi n = l/h\de/dp l x n \ are the longitudinal components of 
the group velocity. [Note that in the conduction band 
(e > 0) and in the valence band (e < 0) the group veloci- 
ties for fixed electron momentum are oppositely directed.] 
The transfer matrix M\ relates the amplitudes ai n and 
bi n on the two sides from the interface according to 



fb ++ \ 






b-+ 1 






\b+Z J 




\ a+- ) 



(9) 



To build this matrix we employ the continuity of the 
electron wavefunction "Ji (0) = (0) and its derivative 

(0) = (0) at the potential step. In this way we find 
that the transfer matrix takes the form 

M x = B- X A, (10) 

with 



A = 



D 
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P++pi 
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In what follows we characterize waves by their corre- 
sponding momentum and amplitude. In the contact, we 
assume that there are right-moving propagating waves 
(p+ + , a++ = 1) which can be transmitted into region II 
or reflected by the potential step Vb back into region I. 
Reflected waves become left-moving propagating waves 
(p~ + , a |_ =^ 0) and evanescent waves decaying to the 



left (p x , a 7^ 0). From this, Eq. ^ becomes 



/b ++ \ 






b-+ 1 


-V, 




\b~+l) 







(11) 



and the transmission coefficient can be found using 



T=l-|a_ + | 2 . 



(12) 
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This definition is the most convenient for the problem 
studied here, since there is only one left-propagating 
mode in region I, whereas there are parametric regimes in 
which two different right-propagating waves exist in re- 
gion II. Below, we discuss how this scheme is applied in 
the four characteristic cases illustrated in Fig. [3j (a) un- 
strained BLG, and (b, c, d) BLG with various angles 
between the principal axis of uniaxial strain and the crys- 
tallographic direction x in Fig. [TJ 

For the unstrained case, using the plane-wave param- 
eters determined in Eq. Eq. pTj ) becomes 




= Mi 



/ 1 

a 

V o 



(13) 



Solving for the wave amplitude a |_ numerically and then 

using Eq. (12), we obtain the transmission probability 
shown in the middle panel of Fig. [3]ja). This reproduces 
the e — > — e asymmetry for transmission of normally in- 
cident electrons (8 = 0), with vanishing T(e, 0) = for 
e < but finite 7"~(e, 0) for e > 0, found in earlier studies 
of BLG j unctions J22H2S This asymmetry can be attributed 
to the different chirality of charge carriers in the conduc- 
tion and valence band. 

For the strain axis oriented at cf> = 0, using the plane 
wave parameters determined in Eq. |6]) and solving for 
the amplitude a |_ in the set of linear equations in (11 1, 



we numerically obtain the transmission result plotted in 
the middle panel of Fig. [3jb) . Our result shows T(e, 8 = 
0) ^ at any |e| < w, as opposed to unstrained case in 
Fig. [3ja) . The difference between transmission at 8 = 
for unstrained and strained BLG can be explained as 
follows. Expanding the Hamiltonian ([!]) in momentum 
around the Dirac points ±p of the strain-induced mini- 
cones [Eq. &\, and keeping only linear terms, we find 
two Hamiltonians valid at |e| <C w, 



±v 



(e*%)t 
e 4( % 



W = 2(Sp x +i6p y ), (14) 



where Sp is a small deviation of the electron momentum 
from ±po, v* = po/(2m) is the effective Dirac velocity, 
and e 1 ^ is a phase factor which determines the position of 
the Dirac points in the momentum plane. By solving the 
Schrodinger equation for each of these Hamiltonians and 
then employing Eq. (11), we can compute the transmis- 



sion probability due to the states in each cone separately. 
The small insert in the middle panel of Fig. [3j^b) shows 
T(0) at e = — 0.2meV for the left and the right Dirac 
mini-cone in red and green, respectively. Transmission 
to the left mini-cone is zero at 8 = and increases away 
from normal incidence, similar to the case of parabolic 
dispersion. Transmission to the right mini-cone, on the 
other hand, exhibits a maximum at 8 = and slowly 
decreases for angles away from 8 = 0, which resembles 
the situation for monolayer graphenepS Therefore, the 
strain-induced mini-cones modify the chirality of the low- 
energy states. 



For the strain axis oriented at <j> = 7r/4, we first verify 
numerically for every angle of incidence and energy which 
momenta correspond to plane waves moving to the right 
(left) and evanescent waves decaying to the right (left), 
respectively. Then we use Eqs. (10) and (11) and solve 

for the amplitude a |_ of the wave reflected back into the 

lead I, taking into account all physically allowed evanes- 
cent and propagating modes in region II. The transmis- 
sion shown in the middle panel of Fig. [3^c) exhibits two 
distinct peaks, as long as the Fermi line is split into two 
pockets. For a small energy range below the LiTr in the 
valence band, T (e,8 = 0) ^ 0, which again can be at- 
tributed to the strain-induced modification of chirality 
of the low-energy states. Beyond the LiTr, where the ef- 
fect of strain becomes weaker and the Fermi line becomes 
circular, we find that T(e, 8 = 0) —> 0. 

For the strain axis oriented at <j> = ir/2 [with plane 
wave parameters determined in Eq. Q], for all transverse 
momenta and energies allowing for propagating states in 
region II the corresponding linear system of equations 
is again the same as in Eq. (13). Solving these equa- 



tions numerically, we obtain the transmission probabil- 
ity shown in Fig. [3jcL) . As a function of incidence an- 
gle 8, the transmission now exhibits two distinct peaks 
for all energies in the considered range. As in the non- 
strained case, Fig. a) , this orientation of the strain de- 
livers T(e < 0, 8 = 0) = 0. 

Irrespective of the modifications of chirality, in all four 
cases there is a marked difference in the transmission 
strength for e > and e < 0. For e < 0, the interface 
which is an n-p junction and an electron incoming in the 
conduction band of the lead (region I) emerges in the 
valence band at the right from interface (region II). For 
e > the electron stays in the conduction band both at 
the left and right from the interface, which is a better 
transmitting n-n' junction. 



B. Conductance 

Based on the above results for the transmission prob- 
ability, we employ the Landauer-Biittiker formalism^ to 
calculate the conductance of the n-p or n-n' junction. 
Taking into account two valleys and two spins, as well as 
integrating over the angle of incidence and electron en- 
ergy (as determined by the Fermi distribution with finite 
temperature T), we arrive at the junction conductance, 



(15) 
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FIG. 3. (Color online) Left: isoenergetic lines at e = —2, —5, —6 and — 8meV for strained bilayer graphene, with w = 5meV. 
Center: transmission probability T(e,8) across a single potential step (n-p or n-n' junction), from a highly doped region to a 
barely doped region, as a function of energy and incidence angle of incoming electrons. Right: linear response conductance of 
the junction as a function of chemical potential /i and temperature T. Results are shown for unstrained bilayer graphene (a), 
as well as uniaxially strained bilayer graphene for various orientations of the strain axis with respect to the crystallographic 
axis x in Fig. [l] = (b), cj> = 7r/4 (c) and cj> = n/2 (d). 



Here, \p ~ 2nhWT/2mVo is the Fermi wavelength in 
lead I, ks is the Boltzmann constant and 7~(e, 9) is the 
transmission at fixed energy and angle of incidence, de- 
termined above. Performing numerical integration we 
obtain the conductance as a function of chemical poten- 
tial and temperature, which is shown in the right panels 
of Fig. § 

For unstrained BLG, Fig. [3]ja), where the dispersion is 
parabolic, the conductance exhibits a minimum at fi = 
and has an asymmetric but monotonic behavior for both 



fi < (n-p junction) and /i > (n-n' junction). Further- 
more, for fi = the conductance increases monotonously 
with temperature. In the strained cases, Fig. c, 
d), the increase in conductance with temperature is still 
seen. However, depending on the strain orientation </>, 
G (/i, T) can be monotonic or non-monotonic. For = 0, 
Fig. [3^b) , the conductance at low temperatures exhibits 
an anomaly — in the region /i < there is an additional 
local minimum, as well as a local maximum. The local 
maximum is located at /i ~ — w (// ss — 5meV for pa- 
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rameters used in the figure), which corresponds to the 
LiTr energy. For = 7r/4, Fig. [3^c) , the conductance 
at low temperatures exhibits a protrusion and a shift in 
slope, which again occur near the LiTr in the valence 
band. For <f> = n/2, Fig. [3^d), the conductance G(fi,T) 
is monotonic and quite similar to that calculated for BLG 
with a parabolic spectrum. 

All conductance plots show an asymmetry about the 
minimum at fi = 0. As discussed for the transmission 
probability, this difference is determined by the chiral 
sublattice structure of the plane waves, which in BLG 
suppresses the transmission at a potential step between 
regions of opposite polarity. Note that the anomalous 
behavior at /i = — w is specific for the n-p junction regime 
of the system, and does not occur in the n-n' junction 
regime, which does not exhibit an anomaly at ji = w. 
These features allow one to single out the anomalous T- 
and /z-dependence of the conductance for junctions with 
different orientation of the strain. 



IV. TRANSPORT ACROSS N-P-N AND N-N'-N 
JUNCTIONS 

All the features found in the parametric dependencies 
of the transmission across a single potential step appear 
also in the transport properties of the two-terminal bal- 
listic device (with two steps) sketched in Fig. [2j In par- 
ticular, the e — > — e, and fi — ¥ —fi asymmetry and the 
anomalous temperature dependence at /x = — w (in the 
vicinity of the saddle point in the valence band) also per- 
sist in this "potential barrier" geometry, and indeed are 
further enhanced. In addition to those, the energy and 
angle dependence of the transmission coefficient acquires 
a resonance structure due to the interference between 
multiply reflected waves (Fabry-Perot resonances). To 
take this into account, we compute the transmission of 
the device sketched in Fig. [2] considering both interfaces, 
as well as the ballistic electron propagation between the 
interfaces. 

The transfer matrix, Mi = B~ 1 A 1 of the first inter- 



face is given by Eq. (10 1. Due to symmetry, the transfer 
matrix of the second interface is Mi = M^ 1 = A~ 1 B. 
The transfer matrix of the whole system (n-p-n or n-n'-n 
junction) is then given by 



S = M 2 SN'h = M^SMi , 



where 
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(16) 



(17) 



describes the ballistic electron propagation inside the 
"barrier" region II. Note that the factors in the matrices 
A and B that normalize the plane wave states in region 
II to normal flux cancel out in the matrix S. From this, 



we can relate the amplitudes a;„ of the wavefunction in 
the source lead, Eq. ([8]), to the amplitudes q„ of the 
wavefunction in the drain lead, 



(x) 




a iPx x 



by 
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To determine the transmission coefficient 

TM) = | C++ | 2 , 

we take boundary conditions c |_ = c = 

a ++ = 1, and find c ++ by solving the equation 
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The numerically evaluated transmission probability T 
is plotted in the left panels of Fig. [4j for the same range of 
angles and energies (— 10meV< e < lOmeV, V a — 50meV 
and w = 5meV) as in Fig. [3} As in Refs. [21 [Ml and [57], 
the presence of two reflective interfaces in a BLG device 
causes the appearance of resonances with high transmis- 
sion. In the figure, these are seen as bright strips. The 
scale of the oscillations becomes finer for a longer sample 
length. 

The right column in Fig. [3] shows the finite- 
temperature conductance in a long sample. The interfer- 
ence fringes are washed out by smearing of the Fermi step 
(at ksT h 2 /(mLz)) and the integration over the an- 
gle. To obtain this finite-temperature conductance, one 
can use the exactly calculated T(e, 8) and insert this into 
Eq. (15). Here we describe an accurate approximation 



of these results, which allows one relate the pronounced 
anomalies of the finite-temperature conductance to angu- 
larly smoothed transmission probabilities (T) (smeared 
over a small angle range SO covering many oscillations), 
shown in the middle column in Fig. [4] Conveniently, in 
the limit of L x — > oo, evanescent modes die off before 
reaching the second interface, so that in region II only 
plane waves (with real p x ) contribute towards transmis- 
sion. To eliminate the negligible contribution of evanes- 
cent waves, we first restrict the analysis to the range of 
angles A8 = 26 c where plane waves exist inside the bar- 
rier; 9 C is energy dependent and different for each orien- 
tation of applied strain. Then, we group the exponents 
which emerge from Eq. (171 into propagating and decay- 
ing waves (where the latter have complex p x ), and for 
decaying waves approximate tanhdlm^JjLj;) — ¥ 1 and 
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FIG. 4. (Color online) Transmission coefficient and conductance of n-p-n and n-n'-n junctions with non-strained bilayer 
graphene (a), as well as strained bilayer graphene with the uniaxial strain axis at an angle for = (b), (f> — n/4 (c) and 
<f> = 7r/2 (d) from crystallographic axis x. Left: transmission probability T(e,9) obtained in an exact calculation. Center: 
transmission probability obtained by averaging over fast oscillations after the contribution of evanescent waves is neglected. 
Right: linear response conductance as a function of chemical potential and temperature. All calculations are performed with 
w = 5meV, Vq = 50meV, and L x — 1/im. 



cosh (2|Im[p x ] \L X ) — > 0. The conductance then follows 
from 



G( M ,T) = 



h \p 
1 1 



.9 



(19) 



de 



-TF*-\e + V \ 



2ir4k B Tj_ oc ^ CQsh 2 
(T(e,e Q ))cos{e )d9 . 



Since the details of the analysis of (T) depend on the 
electron energy and orientation of the strain axis, we 
sketch the derivation separately for the corresponding 
characteristic parametric regimes. 

Firstly, for the range of parameters for which Eq. Q 
has only two real solutions, the described procedure leads 
to an expression of the form 

~ X 2 + X 3 cos (2p++L x ) + X 4 sin (2p++L x ) ' 
Here Xi are non-oscillating functions of p ™, p™, p y , e, 



9 



Vo, w and <p which are not given explicitly due to their 
complexity. To average T(e,6), we first expand the real 
momenta in terms of small deviations 58 in the angle, 
9 = 80 + 58, about some — 8 C < 9q < 8 C , such that 



Px 



08 



As such, 
TMo) 



Xi 



X 2 + X 3 cos ($ + A59) + A 4 sin ($ + A59) ' 



where $ = 2L X p+ + \ g=0o and A = 2L X (0p++ /08)\ e=6o . 
Imposing AA58 = 2ir, the average transmission over one 
period is 



(TMo)) 



1 

2^ 



2tt 



11 X 2 + X 3 cos(z) + A 4 sin(z) 



(20) 



z 



A58. 



Secondly, for the range of parameters where Eq. Q has 
four real solutions, p~ + = —p x + and p x ~ = — p~ _ [such 
as encountered in Fig.^b)], fast oscillations in the trans- 
mission coefficient are due to combinations of sin(p+ n L x ), 
sm(2p+ n L x ), cos(p^ n L x ), and cos(2p+ n L x ). Expanding 
in terms of small deviations in angle, 



v +n \ 



08 



and denoting $„ = L x pt n \e=e and Ai = 
La; (dp x n /d6)\ s _ g , we find that the interference fringes 
are encoded in the factors sin($„ + A n 58) and cos(<I> n + 
A n 58). Inspection of the constant pre- factors reveals 
that A + w A-. Neglecting the phase <fr„ and impos- 
ing A+A58 — 2ir, the averaged transmission over one 
period can then be written as 



and z — A+59; here Xi are non-oscillating functions of 
the same parameters as in the previous cases. The spe- 
cific expressions are again omitted because of their com- 
plexity. 

In the analytical part of the studies of the transmis- 
sion problem, all functions Xi and Xi have been found 
using the symbolic mathematical software. The results 
of the integrals in Eqs. (20l) and (21) are shown in the 
central column in Fig. H The doping (chemical poten- 
tial fi) and temperature dependence of the two-terminal 
conductance of the device follows from Eq. ( 19 ), and coin- 



cides with a high accuracy with the one calculated using 
Eq. (15) together with the exact values T(e,8). 

The behavior of G (/i, T) in the right column of Fig. 4] 
displays all the features of the conductance of a single 
step in enhanced form. In particular, the conductance 
for <ft — 0, Fig. [4jb) , exhibits a local maximum and a 
second local minimum positioned at the same chemical 
potentials as for a single junction. For cf> — tt/2, Fig.BFd), 
the conductance is monotonic. For (f> — 7r/4, Fig. Wc), 
the protrusion in the conductance of a single junction (at 
the LiTr) has developed into a clear local maximum. 



V. SUMMARY 

In this article we have shown that the linear response 
conductance G(/i, T) of an n-p-n junction in strained 
bilayer graphene has a non-monotonic dependence on 
doping and temperature, which varies in size and form 
as a function of the crystallographic orientation of the 
principal strain axis. To understand this behavior we 
studied the transmission and conductance for a single 
interface (n-p junction), and used the obtained results to 
conclude that the non-monotonic behavior is due to the 
change of the chirality (sublattice composition) of the 
electronic plane wave states in the vicinity of the saddle 
point (Lifshitz transition) in the low-energy electron 
spectrum of strained bilayer graphene. 



<TMo)) = 



1 

27 



2jt 



*(*) 

y(z) 



dz. 



(21) 
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